Electronic structure of an electron on the gyroid surface: A helical labyrinth 
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Previously reported formulation for electrons on curved periodic surfaces is used to analyze the 
band structure of an electron bound on the gyroid surface (the only triply-periodic minimal surface 
that has screw axes). We find that an effect of the helical structure appears as the bands multiply 
sticking together on the Brillouin zone boundaries. We elaborate how the band sticking is lifted 
when the helical and inversion symmetries of the structure are degraded. We find from this that 
the symmetries give rise to prominent peaks in the density of states. 



While the physics of electron systems in crystals has 
been firmly established, it should be interesting to con- 
sider electrons bound on infinite, periodic curved sur- 
faces, which we may envisage as crystals composed of sur- 
faces. We have previously calculated the band structure 
of an electron confined on triply periodic (i.e., periodic 
along X, y, z) minimal surfaces (called P- and D-surfaces) 
and found that bands and the Bloch wavefunctions are 
basically determined by the way in which the "tubes" are 
connected into a network. P,Q 

One asset of the crystal of surfaces is that we can 
deform them. We have indeed found that the P- 
and D surfaces, which are mutually Bonnet (conformal)- 
transformed, have related electronic structures. We can 
even twist the tubes, for which there is a special in- 
terest: Schoen pointed out in the late 1960's the fol- 
lowing. P-surface (a cubic network of tubes) and D- 
surface (a diamond network) are typical triply periodic 
minimal surfaces, where minimal surfaces are defined 
as negatively-curved surfaces that have minimized areas 
with the mean curvature (^ (ki -I-K2) with ki , K2 being the 
principal curvatures) vanishing everywhere on the sur- 
face. In addition, however, there exists a third one which 
Schoen called the gyroid(G-) surface. The P-, D-, and G- 
surfaces are related via the Bonnet transformation, where 
the Bonnet angle (a parameter in the transformation) 
/3 = 0,38.0°, 90° correspond to P, G, D, respectively G- 
surface is unique in its triply helical structure, on which 
there are no two-fold axes, nor straight lines, so we may 
call the surface a helical labyrinth. G-surface has been 
explored from the viewpoints of chemistry, crystallogra- 
phy and material science as well. Namely, the gyroid 
structure is known to occur, topologically, in some classes 
of crystal structures, which include clathrate compounds 
(such as Ba6Ge25[3,|l]), zeolite structures (such as MCM- 
and ceramic structures fabricated with copolymer 
templates^. A particular interest, naturally, is how the 
helical structure affects the electronic structure on G- 
surface. 

Here we have calculated the band structure of an elec- 
tron bound on G-surface, adopting the formulation of our 
previous workj^. We find that an effect of the helical ge- 
ometry appears as a set of multiple band sticking phenom- 



ena at the Brillouin zone edges. We have identified that 
the band sticking is related to the existence of screw axes 
and the inversion symmetry by checking that a degrada- 
tion of the helical or inversion symmetries of the structure 
lifts the sticking, in an atomic model ( "graphitic sponge" ) 
realizing the structure. We conclude from this that the 
symmetries give rise to prominent peaks in the density 
of states. 

We start by recapitulating the formulation for repre- 
senting minimal surfaces. We consider a two-dimensional 
surface r{u,v) = {x{u,v),y{u,v), z{u,v)) embedded in 
three spatial dimensions as parameterized by two coordi- 
nates u,v. When the surface is minimal, we can exploit 
the Weierstrass-Enneper representation given as 

r(u,w)=Re( / {I ~ w^)F{w) dw, 

\J W() 

/ i{l + w^)F{w)dw, 2wF{w)dw], (1) 

where w = u -\- iv and F{w) = ie*^L/Vl — Uw^ -\- 
with the Bonnet angle (3 = 38.015° for G-surface, and 
L the linear dimension of the unit cell. As stressed in 
0, Weierstrass-Enneper representation that completely 
defines the surface is specified solely by its poles, which, 
in real space, correspond to 'navels' (umbilical points) in 
differential-geometrical language, at which ki = K2 = 0. 
So a periodic minimal surface has negative curvatures 
everywhere except at the navels. 

Figure n shows a primitive patch of G-surfacc, which 
corresponds to eqn.Q with < 9 < n, — 7r/4 < (j) < n/A 
when we stereographically map {u, v) to a unit sphere 
{0, (j)) with w = u-\-iv = cot(6'/2)e*'^. The full surface is 
depicted in Fig. [21 for its cubic unit cell (containing two 
bcc cells). The surface has 90° helical symmetry axes 
along X, y and z directions, respectively, while the chi- 
rality (right- or left-handedness) is opposite across neigh- 
boring helices. This is how the surface is a network of he- 
lical tubes along x, y, z connected into a single labyrinth. 

While the helical structure is unique to the G-surface, 
it shares with P- and D- surfaces the property called the 
balance surface (or the in-out symmetry), which is defined 
as the curved surface that divides the three-dimensional 
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space into two, congruent labyrinths. For the G-surface 
the divided spaces are mirror inversions of each other. 
The surface itself is not chiral, being symmetric with re- 
spect to flat points (the center of the primitive patch 
in Fig. U) which corresponds to the navels. It has been 
shownp] that a balance surface is characterized by a pair 
of space groups {G,H), where G maps one side of the 
surface to itself or the other side (and one labyrinth to 
itself or the other labyrinth), while H, a subgroup of G, 
maps each side to itself (and each labyrinth to itself). For 
G-surface, {G,H) = (/aSd, /4i32). 

Schrodinger's equation on a curved surface takes dif- 
ferent forms between the following two cases; one is to 
consider electrons bound to a thin, curved slab of thick- 
ness d, where the limit d is taken^, while the other 
is to ignore the perpendicular degree of freedom from the 
outset. We adopt the former as a physical approach, for 
which the equation reads Q 
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where {q^,q'^) = {u,v), gij the metric tensor with sum- 
mations over repeated indices assumed, and ki,K2 the 
local principal curvatures (with ki + K2 = everywhere 
for a minimal surface by definition). An effect of the cur- 
vature of the surface appears as a curvature potential, 
— {h^ /8m){Ki — ^2)^, which has its minima at navels. 

As shown in Ref. 1], the Weierstrass-Enneper repre- 
sentation transforms Schrodinger's equation into a neat 
form of 
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Since the Bonnet transformation preserves the metric 
tensor as well as the Gaussian curvature, the surfaces 
connected by a Bonnet transformation obey the identi- 
cal Schrodinger's equation, as seen in cqn.® where F 
only enters as \F\. However, this only applies to a unit 
patch, so that the band structures are different between 
P,D, and G, since the way in which unit patches are con- 
nected is different. 

So we have solved the equation for the appropriate con- 
nection between the patches to obtain the band structure 
for G-surface as shown in Fig. |21 We can immediately 
notice that the bands are multiply degenerate at some 
symmetric points on the Brillouin zone boundary: six- 
fold degenerate at H-point, four- fold at P, etc, which is 
a phenomenon not seen in P and D-surfaces. The wave 
functions for six-fold states at H are displayed in Fig. 
21 The space group for the G-surface (/aSd) is non- 
symmorphic, i.e., has the screw axes or glide planes, and 
this should cause the degeneracy at the zone boundary, 
generally known as 'the band sticking together' 10] in e.g. 
crystals of Se and Te with helical atomic arrangements. 
We shall elaborate this point below. 




FIG. 1: A primitive patch for G-surface. Its center corre- 
sponds to the navel (see text), around which the surface is 
point-symmetric. 




FIG. 2: A cubic unit cell of G-surface in birds-eye(a) and 
top(b) views. 



We can note in passingthe following: We have pointed 
out in the previous work[lj that Bonnet-connected P- and 
D-surfaces, despite the different band structures, have a 
common set of band energies at special fc-points (Bril- 
louin zone corners, edges, or face-centers). This occurs 
when a wave function on a unit patch can be continued 
as a "tiling" of patches on both surfaces, where the sim- 
plest case is the ground state at F point. In G-surface, 
this occurs at (and only at; a peculiarity of the spiral 
structure) F at which the band energies accommodate 
those of P-surface at F and H points and D-surface at F 
and R as displayed in Fig [3 

A question from the materials science point of view is: 
can we realize the G-surface in some materials? As men- 
tioned above, there are a few classes of materials that 
possess the same space group as the G-surface. One is a 
class of clathrate compounds of group- IV elements (e.g., 
Ge), where Ge2o clusters that include another element 
(e.g., Ba) are stacked in a triply-helical fashion. A con- 
ceptually simpler system would be a three-dimensional 
labyrinth of graphite sheets that forms a triply peri- 
odic surface. Fujita and coworkers have considered such 
systems, and called them "graphitic sponges" Til . The 
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FIG. 3: Band structure of the G-surface displayed on the 
Brillouin zone for the bcc unit cell. The numbers in the paren- 
theses represent the degeneracy of the band sticking. 




FIG. 4: The wave functions in the unit cell of G-surface for 
the six-fold states at the lowest level in H-point (in Fig. 
with positive (negative) amplitudes color coded in red (blue) . 
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FIG. 5: Coincidence of the band energies among G-, P- and 
D-surfaces. 



(a) Ia3d 



(b) Ia3 (c) I4i32 





-2,8 





i'll 








.■!»'! 






E 



















(3) 






(3) 












fi) 


c 










ii! 






(1) 




; b; 


(1) 




-■:;;;A 


(3) 






..■■■'a 


(3) 




FIG. 6: (a) Low-energy band structure for the graphitic 
G-sponge, to be compared with Fig. |21 constructed from a 
graphite fragment shown at the top. We take here the hop- 
ping parameter t as the unit of energy and the tight-binding 
band center as the origin, (b) The band structure when we 
introduce a potential that breaks the helical symmetry, or (c) 
the inversion symmetry. Extra potential introduced to de- 
grade the symmetry is shown at the top, where an open circle 
represents the potential for every other patch, while filled one 
for all the patches. The panels (d) and (e) show the density 
of states for the systems (a) and (b), respectively, in units of 
1/iVt) {V: unit ceU volume). 
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negative-curvature fuUerene (or Cgo zeolite) [12] is indeed 
a realization of the P-surface if we smear out atoms into a 
surface, [l^j The band structure of atomic networks such 
as the Cgo-zeolite is expected to basically reflect the prop- 
erties of an electron on the curved surface, as far as the 
effective-mass formalism is applicable and effects of the 
odd-membered atomic rings are neglected. This has in 
fact been shown for P- and D-surfacesQ- Fabrication 
of graphitic sponges has been experimentally attempted 
with a zeolite as a template so we should end up 
with a G-surface sponge if we use a G-structured zeolite 
(MCM) as the template. 

The electronic structures of the negative curvature car- 
bon networks, including G-surface-like structures, were 
investigated with a first principles method in Ref.'lS'l. 
However, some of the G-surface-like models (E and O 
in that paper) contain seven- membered rings, which de- 
grade the symmetry below that of the G-surface, and 
the band degeneracy was not discussed, either. Here we 
examine the band structure of the graphite sponge that 
shares the global topology and the symmetry with the G- 
surface to investigate how the degeneracy would be lifted 
when the symmetry is degraded. We constructed the net- 
work by arraying the unit patches displayed in Fig. IBfa) 
(corresponding to Fig. where the resultant structure 
contains no odd-membered rings, and the band structure 
is calculated in the one-band, tight-binding model with 
the hopping parameter t. 

The result in Fig. Elb) shows that the bands in the 
low-energy region have a one-to-one correspondence with 
those for the G-surface, as expected, so the electronic 
structure there is basically determined by the structure 
of the surface on which the atoms reside. Then we have 
degraded each of the two symmetries (helical and inver- 
sion) , by introducing an extra potential. We first destroy 
the helical symmetry in Fig|H^b), with the space group 
laSd reducing to /a3, with the extra potential introduced 
on every other patch. We next destroy the inversion in 



FigEIc) (with the space group reducing to /4i32) with 
the extra potential on triangular sites around the center 
of every patch. The amplitude of the extra potential en- 
ergy is taken to be 0.3t or 0.5i in (b), (c), respectively, 
which are small enough for the correspondence among 
the bands to be retained. The band structures show that 
the multifold stickings split into lower degeneracies, al- 
though some degeneracies in (c) (with the inversion sym- 
metry broken) are lifted only slightly. This endorses that 
both the helical and inversion symmetries are essential 
for the band sticking. The panels (d) and (e) display 
the density of states (DOS) before(d) and after(e) the 
helical symmetry is broken. We recognize that the first 
and second prominent peaks (labeled as A and B, cor- 
responding respectively to six-fold and two-fold stickings 
in (a)) split with significantly reduced peak heights after 
the symmetry is degraded. The effect is more noticeable 
for the peak B, which is due to the splitting of nearly 
flat bands. The density of states for the broken inversion 
symmetry (not shown) exhibits a similar behavior. So 
the band sticking due to the G-surface space group can 
generally give rise to prominent peaks in the density of 
states. The third peak (C) is also suppressed, but this is 
caused by a change in the band curvature. 

Let us finally mention a possible relevance to real 
materials. A recent band calculation for the clathrate 
Ba6Ge25, which shares the space group with G-surface, 
shows that the DOS peak at the Fermi level becomes 
split when the structural symmetry is reduced, which 
they suggest causes a band Jahn- Teller instability, and 
should be responsible for an experimentally observed 
phase transition^^. Our result suggests that materials 
with the G-surface symmetry generally have a potential 
for such instabilities and corresponding phase transitions. 

We wish to thank Shoji Yamanaka and Yasuo Nozue 
for illuminating discussions in the early stage of the 
study. 
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